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We present explicit wavefunctions for quasi-hole excitations over a variety of non- 
abelian quantum Hall states: the Read-Rezayi states with A: > 3 clustering properties 
and a paired spin-singlet quantum Hall state. Quasi-holes over these states consti- 
tute a topological quantum register, which can be addressed by braiding quasi-holes. 
We obtain the braid properties by direct inspection of the quasi-hole wavefunctions. 
We establish that the braid properties for the paired spin-singlet state are those of 
'Fibonacci anyons', and thus suitable for universal quantum computation. Our deriva- 
tions in this paper rely on explicit computations in the parafermionic Conformal Field 
Theories that underly these particular quantum Hall states. 

1 Introduction 

The realization that quantum Hall systems may harness what is called non-abelian braid 
statistics has led to two exciting prospects. The first is that experiments can be set up where, 
for the first time, the existence of non-abelian statistics in nature can be established. The 
second prospect, following early ideas of Kitaev [Tj, is that systems exhibiting non-abelian 
braid statistics can give rise to quantum registers and may offer unique possibilities for 
what has come to be known as topological quantum computation or fault tolerant quantum 
computation. 

Proposals for experimental detection of non-abelian statistics |2l El IH have mostly 
focused on the so-called = 5/2 quantum Hall state, which is believed to be described 
by the 'pfaffian' or Moore- Read state [HlITj. To this date, this is the single quantum Hall 
plateau where the indications for an underlying non-abelian state are rather firm. For the 
plateau at v = 12/5 there is a more speculative case for a connection with a non-abelian 
k = 3 Read-Rezayi state jH]. Proposals to test if the state observed aX u = 12/5 is indeed 
non-abelian have been described in P| llOj. 

An interesting domain for the search for non-abelian quantum Hall states is that of 
multi-component quantum Hall states, which can be realized as double layer or as spin- 
singlet states. The double layer case in particular may offer the experimental flexibility 
needed for tuning into a regime where a non-abelian state takes the upper hand. Examples 
of candidates are the spin-singlet (or double layer) states of [TI], exhibiting a separation of 
spin and charge, and the spin-singlet analogues of the Moore-Read and Read-Rezayi states 
that we introduced in jT2] and studied further (with N. Read and E. Rezayi) in |13J. 
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Ideas for topological quantum computation in quantum Hall systems boil down to the 
following. Having available a non-abelian quantum Hall state as the underlying medium, 
the injection of quasi-hole excitations is known to open up an internal space (the 'quantum 
register'), whose dimensionality grows with the number of quasi-holes. This quantum register 
can be addressed by performing adiabatic quasi-hole braidings, which give rise to matrices 
acting on the register. 

It is well-known that the non-abelian braiding in the Moore-Read state is not suffi- 
ciently rich to enable universal quantum gates on the quantum register (note that there are 
proposals to combine topological and non-topological operations to obtain universal gates 
in this state ^ll^^j). It is also known that 'higher' non-abelian quantum Hall states, 
such as those in the Read-Rezayi series with A; = 3, 5, 6, . . ., do offer a prospect of universal 
quantum computation [TTl HH] , in the sense that universal quantum gates, such as the 2- 
qubit CNOT gate, can be approximated with arbitrary precision by a well-chosen sequence 
of braid matrices fnil2ni- In this paper, we shall establish that the paired spin-singlet state 
(AS state) of [12] is also universal for quantum computation. It is the unique paired quantum 
Hall state with this property. 

The simplest scenario for braid matrices with sufficient structure for universal topolog- 
ical quantum computation is offered by the so-called 'Fibonacci anyons'. In abstracto, these 
are particles of two types, '0' and '1', with fusion rules 

0x0 = 0, 0x1 = 1, 1x1 = + 1. (1.1) 

Through a universal connection between, on the one hand, fusion rules and, on the other, 
braid matrices, it has been established that the simplest braid matrices for Fibonacci anyons 
are, in the notation of 



R 



(-l)"3/5 



T —T 



B =(-!)-"'{ (^i,:,/.^ ) . 



r 



where r = \{V5 — 1). 

Quasi-holes over non-abelian quantum Hall states cannot straightforwardly be iden- 
tified with 'non-abelian anyons', but their are important parallels, in particular where the 
fundamental relations between fusion and braiding properties are concerned. These relations 
were first studied in the context of the algebraic approach rational conformal field theories, 
as developed by Moore and Seiberg in |22i|. In the context of quantum Hall systems, these 
same relations have been exemplified in the work of Slingerland and Bais who used 
an associated quantum group structure to obtain explicit results for braid matrices for the 
Read-Rezayi quantum Hall states for general k. 

Before the work of [23j, Nayak and Wilczek [21] had derived explicit wavefunctions 
for four quasi-hole excitations over the Moore-Read state. These wavefunctions provide de- 
tailed information on the internal state associated with four quasi-holes, as a function of 
the locations of these excitations. From these wavefunctions, braid properties are derived by 
direct inspection. The work of [21] was based on a 'coordinate' rather than an 'algebraic' ap- 
proach to (rational) Conformal Field Theory (CFT), and employed bosonization techniques 
for mastering the relevant CFT correlators. 
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It is the purpose of this paper to present exphcit expressions for wavefunctions of 
quantum registers associated to non-abehan quantum Hall states that are sufficiently rich 
to enable universal topological quantum computation. We will in particular focus on two 
distinct quantum Hall states that both give rise to braid matrices of the type displayed in 
eq. ()1.2|1 (up to additional abelian phase factors). The first is the k = 3 Read-Rezayi state 
and the second is the paired {k = 2) AS spin-singlet quantum Hall state. We shall also 
write some of the wavefunctions for quasi-holes over the general k Read-Rezayi states. Note 
that throughout this paper, we will assume that the only effect of braiding comes from the 
explicit monodromy. 

The appearance of the 'Fibonacci-type' braid matrices in the quantum Hall systems 
can be understood from a coarse-graining of the fusion rules of the parafermionic CFTs 
underlying these states. For the example of the spin-singlet state, this takes the form 

'0' = {l,Vi,^2,V'3} , '1' = {^1,^2, as, p} (1.3) 

with the ipi denoting the parafermion sectors and cTj, p labeling the various parafermion 
spin fields in the CFT. The fusion rules of these fields (see table E] in section IB.ljl are such 
that the coarse graining into '0' and '1' respects the relations given in eq. (jl.lll . We would 
like to stress that the 'Fibonacci anyon' aspect of the quantum Hall quasi-holes captures a 
limited fraction of their relevant properties: for example, the fundamental quasi-holes come 
with different sets of quantum numbers and the detailed fusion rules (and Operator Product 
Expansions) of the fields in the parafermionic CFTs lead to detailed structure expressed in 
the wavefunctions that we derive in this paper. 

The 'coordinate CFT' approach that we follow to derive the quasi-hole wavefunctions 
delves deep into the results for Wess-Zumino-Witten (WZW) and parafermion CFTs derived 
in the mid-1980s. We shall in particular rely on results of Knizhnik and Zamolodchikov (KZ) 
(2n] • Their expression for 4-point functions of particular primary fields in WZW models will 
form the cornerstone of the 'contraction arguments' that we employ to determine closed form 
expressions for the various quantum register states that we consider. The 'master formulas' 
that we develop enable an easy evaluation of correlators that have until now not appeared in 
the literature, and that are not easily computed with the methods of KZ. We present some 
explicit examples in sections I A. 41 IB.41 

The presentation in this paper is organized as follows. In section |21 we will explain 
the method we use to obtain the quasi-hole correlators by using the Moore-Read state as 
an example. In sectional we apply this method to the k = 3 Read-Rezayi states. We will 
provide a detailed derivation of the quasi-hole wavefunctions, and use these to calculate the 
braid behaviour of the quasi-holes. In section |3J we apply our method to the paired spin- 
singlet states proposed by the authors. In sectional we compare the braid results of sections 
El and m with the results obtained by using the associated quantum groups. In the two 
appendices 1X1 and iBl we give the details of the parafermion CFT's corresponding to su{2)k 
and su{3)2 respectively. Namely, we provide the fusion rules, the details of the Operator 
Product Expansions, including the various coefficients, the spin-field correlators used to 
derive the quasi-hole wavefunctions and the braid behaviour of the various hypergeometric 
functions which show up in the correlators. In addition, we give some new parafermionic 
correlators, which were used to obtain the various OPE coefficients. In section IA.61 we give 
the quasi-hole wavefunctions for the Read-Rezayi states for arbitrary k. 
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2 General form of quasi-hole wavefunctions 

In our analysis of non-abelian quantum Hall states, we rely on the so-called qH-CFT con- 
nection where wavefunctions of a quantum Hall (qH) system are expressed as chiral correla- 
tors (conformal blocks) of an associated Conformal Field Theory (CFT) The connection 
hinges on the identification in the CFT of an electron operator ipdz) , carrying charge q = — 1. 
The quantum Hall wavefunction is then expressed as 

^'qH(^l, . . . ,^Ar) = lim {Z^)~ {tp^{zi)iJe{z2) ■ ■ ■i^c{zN)Qhg{Zoo)) (2.1) 

Zoo CO 

with Qbg representing a neutralizing background (ionic) charge and the factor (-2oo)~ is 
included to obtain a non- vanishing result. (The case of spin-full fermions has some additional 
structure.) Note that we drop the Gaussian factors throughout this paper. 

The injection of a quasi-hole at position w is represented by the insertion in the CFT 
correlator of the quasi-hole operator 0qh(w). Mutual locality of the quasi-holes and the 
electrons in the quantum Hall condensate implies that the Operator Product Expansion 
(OPE) between electron and quasi-hole operators is of the form 

UMM^) = {w- zr%^{z) (2.2) 

with n a non-negative integer. Note that the mutual locality puts a constraint on the possible 
quasi-hole operators. 

In all cases studied in this paper, the electron and quasi-hole operators are expressed 
in terms of free bosonic fields (representing charge and spin) and of a parafermionic CFT. 
The latter are closely related to Wess-Zumino-Witten (WZW) theories. For the order-fc 
Read-Rezayi (RR) states the 'CFT-data' are: SU{2)k WZW theory and the associated 
parafermions, while the paired spin-singlet states are connected to the SU{?))2 WZW theory 
and to the associated higher-rank parafermions. We refer to [2^] for a review and further 
details. 

For the Moore- Read (MR) state (which is the k = 2 member of the RR series), the 
parafermion theory reduces to a single real (Majorana) fermion ip{z). This allows a di- 
rect evaluation of the electron wavefunction using the Wick theorem, leading to a 'pfaffian' 
wavefunction. Quasi-holes over this pfaffian state can be characterized by pair braking in 
the pfaffian BCS factor [27|. The full dependence of a multi-quasi-hole wavefunction on all 
coordinates is set by a CFT correlator 

{a{w,) ■ ■ ■ a{wn)^{z,) ■ ■ ■ ^{zN)f''^ (2.3) 

with a{w) the spin-fields for the Majorana fermion. In a 1996 paper, Nayak and Wilczek 
used bosonization techniques to derive an explicit expression for the full four quasi-hole 
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wavefunction [21]. It has the general form 

A^'^'\{w})^uM^}, {4) + B^'''\{w})^,,M^}, {z}) . (2.4) 

In this expression, the factors \l/i2,34 and \l/i3,24, which are polynomial in all coordinates 
({if}, {2;}), represent independent (in this case: two) ways in which four quasi-holes can 
break up pairs in the electron condensate. The prefactors A^^'^^ and B^^'^^ are given by (note 
that we are giving the result for the fermionic case M = 1) 

A^'\{w^}) =t^(^i2ti^34)^a;^ (^(-l)yi-v^+^l + v^^ 

i?(P)({«;.}) = - iz^(«;^2«^34)ia;-|(l _ x)h (^-^1- + (-1)^1 + , (2.5) 

where p = 0,1 and we used the anharmonic ratio x = ^^^^^^ , (where Wij = {wi — wj)), which 
is the same as the anharmonic ratio used in but differs from the one used in The 
labels (0, 1) refer to the fusion channel of the 4-quasi-hole state; it is this index which acts 
as the qubit index in the context of topological quantum computation. 

The formula p.4|l admits an elegant interpretation: it describes precisely how the fusion 
channel basis of the 4-quasi-hole states decomposes over a basis set by patterns in quasi- 
hole induced pair-breaking in the condensate. This decomposition is given as a function 
of the quasi-hole coordinates Wi, meaning that it can be followed as quasi-particles move 
through the condensate. This in particular implies that quasi-hole braiding properties can 
be read off from these wavefunctions, as was done in Clearly, the information stored 
in these wavefunctions goes well beyond braiding properties; we expect that some of this 
additional structure will be relevant for the optimal design of experimental protocols aimed 
at demonstrating non-abelian statistics and at quantum computation. 

In this paper we show that quasi-hole wavefunction for non-abelian quantum Hall 
states with potential for universal topological quantum computation can be cast in a form 
similar to ()2.4j) . To achieve this goal, we rely, in a first step, on known expressions for the 
multi-parafermion correlators representing the quantum Hall condensate in the absence of 
excitations (|H1 EHl ^1 EH])- Analyzing the factors associated with the injection of quasi- 
holes (of various kinds) then leads to 'master formulas' not unlike ()2.4|) . In a final step we 
consider various 'contractions' of this master formula and use those to relate the coefficients 
such as A^^'^^ and B^'^'^^ to correlators having just four parafermion spin-fields. The latter can 
extracted from j2Sl, where they have been expressed in terms of hypergeometric functions. 

3 Quasi-hole wavefunctions for the k = 3 Read-Rezayi 
states 

In this section, we use the approach outlined in the previous section to obtain the wave- 
functions of quasi-holes over the k = 3 Read-Rezayi (RR) states, jH], which can be viewed 
as clustered analogues of the paired Moore- Read (MR) state [S^. The RR states have been 
studied in detail by various authors. Among the advances are explicit wavefunctions in terms 
of the electron coordinates, both with and without the presence of quasiholes. 
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Even though the structure of the quasi-hole wavefunctions in terms of the electron 
coordinates is known, see [2H1 EOl; the full quasi-hole wave functions, which also exhibit 
the full dependence on the quasi-holes coordinates, have only been known for up to four 
quasi- holes in the MR state (see |^) and for an arbitrary number of quasi-holes in the 
Laughlin states. In this paper, we fill in this gap, and calculate the full (four) quasi-hole 
wavefunctions for the RR states and their spin-singlet analogues. In this section, we provide 
the details of the k = 3, M = case of the RR states. We give the general k, M results 
in the appendix ^ where we also give the details of the fusion rules, OPE's, some general 
parafermion correlators, and the braid relations. 

Before we turn to the quasi-hole wavefunctions, we first give the wavefunction of the 
k = 3 RR states without quasi-holes. In this case the number of electrons (even though for 
M = 0, the particles are bosons, we will refer to them as electrons) has to be a multiple of 
3. It was shown in [2H] that the following wavefunction is equivalent to the wavefunction 
presented in the original paper [Sj. Divide the electrons into three groups Sa, a = 1,2,3 of 
equal size. For each group, we write a Laughlin factor 

^L({4)= n (3-1) 

i<j 

To obtain the RR wavefunction, we sum over all inequivalent ways to divide the electrons 
into three groups of equal size 



^'kkHU}) = ^ E [^^.^I^^Is] ' (3-2 



AT 

Si, 82,83 

The normalization is A/" = 3"^/3!, chosen consistently with the operator product expansion 
of the parafermion fields. In effect, the sum amounts to symmetrization of all coordinates. 

The k = 3 clustering property is manifest from the equation (j3.2p : from the wavefunc- 
tion ()3.2|) it is clear that we can put three electrons at the same location, without obtaining 
zero, because there will always be a term in the summation for which the three electrons 
belong to different groups. Putting four or more electrons at the same location gives a 
vanishing wavefunction. 

3.1 The CFT formulation 

The k = 3 RR wavefunctions for electrons and n quasi-holes can be expressed in terms of 
a parafermionic correlator in the following way 

(o-i(Wi) ■ ■ ■(Ti(w„)?/'i(zi) ■ ■ ■ Vl(2;7v))^°'^^ 

X n^-^* - ^jf'^^'Yliw^ - Yliw^ - W,)-^-^^ . (3.3) 

i<j i,j i<j 

From the fusion rules, it follows that the number of electrons and n have to satisfy 
the relation 2N + n = mod 3 in order for the correlator to be non-zero. In addition, a state 
with only one quasi-hole is impossible. 
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3.2 The quasi-hole wavefunctions 

We focus on the case of four quasi-holes. In this case, there are two fusion channels (labeled 
by (0) and (1)) for the parafermion correlator, namely (Xicri ~ ipi fo^' the (0) channel and 
ciai ~ (72 for the (1) channel. 

Following ^.28.^ and specifying to = 3r + 1 electrons, with r an integer, we can write 
the following ansatz for the wavefunction for four quasi-holes, ()3.3p . 

A^'^'\{w})^uM^}, U}) + B^'''\{w})^r,M{^}, M) ■ (3-4) 

Throughout the paper, we will choose the phase of wavefunctions in such a way that the 
function A^^\{w}) has no phase. 

To specify the functions \E'i2,34 and \E'i3,24, we divide the electrons in three groups in 
such a way that contains (A^ — l)/3 + 1 electrons and 5*2 and S3 contain (A^ — l)/3 
electrons. 

Splitting the four quasiholes into two groups, we have 



^13,24 = ^ Yl [n^^^ - ^l)^^^ - n - ^2)(^, - W^)^l^^l^^%'^ , (3.5) 

Si,S2,S3 ie52 jeSa 

with A/ = 3^. In the following, we will use the case of = 4 electrons to determine the 
functions A^^'^\{w}) and B^^'^\{w}). Note that in the functions ()3.5|1 . the quasiholes do 
not have a zero with electrons of the first group before symmetrization. Nevertheless, the 
maximum degree of each z is the same, because the number of electrons in the first group is 
one bigger in comparison to the other two groups. 

There is a third possible way of dividing the four quasiholes into two groups. However, 
this does not give an independent function, because we have the following relation 

^14,23 = X ^12,34 + (1 - a:)^i3,24 • (3.6) 



Relations of this kind were first studied by Nayak and Wilczek in p^, in order to reduce 
the overcomplete basis to a linear independent one. 

We should note that an (explicitly) independent basis for the Read-Rezayi states (for 
arbitrary k) was recently formulated by Read [20] (building on results presented in t32j). 
However, for our present purposes, namely deriving the full quasi-hole wavefunctions and 
studying the braid behaviour of the quasi-holes, it is more convenient to use the basis states 
(|3.5p . because the formulas and the transformation properties under the braiding of quasi- 
holes are simpler. 

The strategy to obtain the functions A'-^-' and -B*-^-* in ()3.4p is as follows. Making use 
of OPE's, including the OPE coefficients, we reduce the correlator in the wavefunction ()3.3|) 
to correlators involving just four o"i^2 fields. The latter can be extracted from the results in 
|25j . In particular we consider the following limits 

(/) Z2 Zi, Za Z3, Z-i Wa, Zi W2, 
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(//) Z2 Zi, Z4 2:3, Z3 W4, Zi ^W3. (3.7) 

In both limits, we first fuse tlie two pairs of tpi fields to two ^/'2 fields. Fusing these 1^2 fields 
with (Ti fields, we obtain (T2 fields (see section \Ka] for more details about the fusion rules). 
Thus, we obtain the following correlators 

lim ((Ti(Ti(7i(Ji'l/'i'l/'l?/'l?/'l)^°'"^^ oc (o-iO"20"iO"2)^°'"^^ 

h.m{aiaiaiai^lJi^lJiipiipi)''^'^^ oc (criai(T20-2)^°'"^^ (3.8) 

Keeping track of the various OPE coefficients and various phases is crucial in this 
procedure. Many of these coefficients for the parafermion CFT are given in [HI] . The other 
coefficients we need were obtained from consistency relations on the four point correlators. 
We give these coefficients in appendix IA.21 

With respect to the various phases which need to be taken into account we note the 
following. To be able to fuse the fields 1^2 with the appropriate spin field, one has to move 
the field ip2{zi) in between the cxi fields. The resulting phase depends on the fusion channel, 
and can be obtained from the OPE of 1JJ2 and the fusion channel under consideration. As an 
example, in the limits above we need to find the phase associated with (Ji{w3)ai{w4)ip2{zi) = 
{—I)°'ip2{zi)ai{w3)ai{w4). In the (0) channel, we have cticti ~ ipi, and the phase a = — | 
follows from 

while in the (1) channel we have (Jicti ~ (T2 and the phase a = — | follows from 

a2{z)ij2{z') ~ e{z'){z - z')-'' ~ (-l)-^£(/)(/ - z)"^ ~ (-l)-^2(^')^2(^) . 

Taking the phases and OPE coefficients into account, we obtain the limits (/) and (//) 
of the wavefunction ^Prr, with p = 0, 1, eq. ()3.3j) 

lim^'J[^ = -(-l)^^(^i'l2W34)^(l - X)^ {(Tl{Wi)a2{w2)(Ti{w3)a2{Wi))^^^ {WI2WUW32) 

(I) y 

lini^'i^^ = {-lY^{Wi2W34)^X~^l -X)t {cri{wi)(Ti{w2)(T2{w3)(T2{w4))^^^ (^34^^14^32) , 

(3.9) 

where we used the notation Wij = {wi — wj) and the following convention for the anharmonic 
ratio X 

^ ^ (W12)(W34) l-x= (^13)(W42) X ^ {Wi2){wu) 

(W14)(W32) (W14)(W32) 1-X (Wl3)(w42)' 

which is the same convention as used in [25J, but differs from the one used in jSIj. The 
explicit form of correlators (cri(j2(Ticr2)^°'^^ and {(7iaia20'2)^^'^^ can be extracted from the 
results of [23]. In appendix IA.3I we present formulas expressing these correlators in terms of 
hypergeometric functions. 
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Because we can easily take the limits (/) and (//) of the functions ^1^12,34 and \E'i3,24 in 
the case of four electrons, namely 

4 

lim ^12,34 = -77^14^32(^42)"^ lim ^12,34 = 

(I) 9 (//) 

4 

lim ^13,24 = lim ^13,24 = -^10141032(1034^ > (3-11) 

(I) (11) '9 

we finally obtain relations for the functions A*^'^'^-' and i?^^'^) which are easily solved to give 

= {wi2W34)^^X^^T^'\x) 

= -(wi2W34)^a;-i^(l - x)T^^\x) 
A« = _(_i)fc7(w;i2w;34)^a;^J-J'^(x) 

= {~1)Ic{wi2W3a)^x-^{1 - x)J^^^\x) . (3.12) 



Here, C = \J r(^)r3(t) ^^"^ functions , with p = 0, 1 and z = 1, 2 are given in ()A.12|) . 

With the explicit form of the functions A^"* and B^^\ we have completely specified the 
wavefunction for the = 3, M = Read-Rezayi state in the presence of four quasi-holes. 
We will give the general fc, M results in the appendix ^ As an immediate application, we 
can study the behaviour under exchange of quasi-holes, by making use of transformation 
properties of the hypergeometric functions, which are also given in appendix 1X1 

3.3 Braid behaviour 

To study the braid behaviour under the exchange of quasi-holes, we first note that the 
anharmonic ratio transforms as x ^— > for (1 <-> 2), a; h-> 1 — x for (2 ^ 3) and a; 1— > ^ for 
(1 3). In addition, we find that, by making use of ()3.6|) . 

^12,34 ^ X ^12,34 + (1 - a;)^l3,24 

^^13,24 ^a;^i2,34 + (l-x)^ 13,24 5 (3.13) 

while the other transformations of the \E''s are clear. The braid transformations of the 
functions JF.*-^-* are given in appendix IA.51 Combining all the results, we find the following 
braid behaviour for general M 



(-1) 



3 



< — (Ci3)?< Ur.^i-lf'-^ij^ 7_^J'iV) (3.14) 



1 
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where we use r to denote the inverse of the golden ratio, i.e. r = ^ ^ . These matrices are 
unitary, and satisfy f/13 = f/23 ■ Uvi ■ U23 ■ Note that if we evenly distribute the phases of 
the off-diagonal elements, which amounts to a 'gauge' transformation, the matrices Uu, U23 
and Ui3 become proportional to the (inverses of the) R, F and B matrices of the 'Fibonacci' 
anyons, as displayed in eq. ()1.2|) . 

4 Quasi-hole wavefunctions for a paired spin-singlet 
state 

4.1 A paired spin-singlet quantum Hall state 

In the search for a topological quantum liquid suited for universal quantum computation, 
the paired spin-singlet quantum Hall state of [12] imposes itself as a natural candidate. This 
state is the k = 2 member of a series of spin-singlet states introduced and studied in [T^ITH]. 
In many ways, these states are direct extensions of the Read-Reazayi states to spin-full 
(spin-1/2) fermions. 

The simplest fermionic spin-singlet state (with M = 1) has filling fraction u = 4/7. 
At this particular filling spin-singlet quantum Hall states have been observed ||M|, but their 
precise nature has never been determined. 

The braiding properties of the paired spin-singlet state are essentially more complicated 
than those of the paired spin-polarized (Moore- Read) state. In fact, we will show that the 
braiding in the paired {k = 2) spin-singlet state is similar to that in the k = 3 Read- 
Rezayi state, the similarity being due to what is known as level-rank duality between the 
affine Lie algebras 5^(2)3 and 5^(8)2. With this, the k = 2 state offers the perspective 
of universal topological quantum computation in a paired quantum Hall state. This being 
enough excitement for us now, we shall in this paper not address the k > 2 spin-singlet 
states. Their wavefunctions can be obtained using the methods described in this paper. 

4.2 The CFT formulation 

The various wavefunctions for the k = 2 spin-singlet state in the presence of quasi-holes are all 
expressed as correlators in a CFT. For M = the CFT is precisely the (chiral) SU (8)2 WZW 
model while for M 7^ we have a deformation thereof, with a modified compactification 
radius of the charge boson yj^, see ^j. In all cases, the theory is conveniently represented 
as a product of the 5^(8)2 parafermions, as introduced by Gepner in [SI] and the CFT of 
spin and charge bosons and ipc- 

The fundamental quasi-holes over this quantum Hall state come in different types. One 
type has spin-1/2 and charge l/(4M-|-8); a second option is to have spin-less quasi-holes of 
charge 2/(4M + 8). For M > the latter have the smaller scaling dimension and are thereby 
the most 'relevant' in the sense of scaling arguments. We expect that experimental protocols 
for the detection of non-abelian statistics and for quantum computation will be most easily 
implemented using the spin-less quasi-holes. 

Using the quantum Hall - CFT connection, we can write the wavefunction for a state 
with the number of quasi-holes of the various types specified as n^, nj^, n^. These numbers 
satisfy 

N^ + n^ = N^ + ni , (4.1) 
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in order for the state to be a spin-singlet and 



3(A^t + A^l) + (n| + + 2^3 = 0mod4 (4.2) 

such that the fields can be fused to the identity sector. The case of only one quasi-hole is an 
exception, because a state with only one quasi-hole is impossible. 
In full generality the quasi-hole wavefunction reads 



*As(^I> ■■■,wl^^wi,...,wi^]Wi,...,Wn:,-,zl,...,zl^-,zi,...,zl^^) = 

hJ hj ij h3 

n-t- 1 M -| — |- I 1 M -| — |- I I 1 M 

{Wl-W!j)'i 3C4Af+3) 11(^1 -wj) 6 3{4M+3) _ y;j)3 3(4M+3) 

i<j i,j i<j 

nt 1 2M -p-|- I 1 2M -p-|- 1 4M 

(W. -Wj)'' 3{4M+3) _ y;^.)6 3{4M+3) | | (y; . _ y;^. ) 3 3{4M+3). (4.3) 



X 



Below we present explicit formulas for these wavefunctions in the special cases of four quasi- 
holes with (i) = 4, (ii) = = 2, (iii) = 4 and (iv) = 3 , = 1. To avoid clutter, 
we write the expressions for M = 0; the braid relations will be specified for general M. 

4.3 Evaluating quasi-hole wavefunctions 

We first give the wavefunction without any quasi-holes in the form given in j^H], see also 
|28j . Assuming = both even we have (the sum is over all independent ways of 
dividing the electrons in two groups, both containing spin-up electrons and iVj/2 spin- 
down electrons; in effect this amounts to symmetrization over the spin-up and spin-down 
electrons) 

*A>r° = ^E*"'(^.'.4)*S'(4.4). (4.4) 



with J\f = 2 2— /2 and with "^f^ denoting the 221 state restricted to zl, Zj, G Sa 

m\^i4) = n - n (4 - 4)' n (4 - 4) • (^-s) 

i<j i'<j' 

iJ'^Sa i',j'eSa i,j'&Sa 

Note that from now on, we will put a prime on the index of spin-down particles. In some 
cases we will drop the arrows on the quasi-hole coordinates. 

The validity of expression ()4.4j) is most easily understood from the characterization of 
the paired spin-singlet as the maximal-degree zero-energy eigenfunction of a specific 3-body 
hamiltonian. 
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4.3.1 The case ^3 = 4 

The wavefunction takes the form 

^As [S333]{wi,W2,Ws,w^; zl, . . . , zj^^; z[,, . . . , zj^, ) = 

{a3{wi)a3{w2)a3{ws)as{w4)'ilJi{zl) . . . V'i(^jv^)V'2(4) • • • V'2(^Jv/)) 



X 



1/2 



^1 ~wjpYli^i'-^jp 



i<3 



(4.6) 



Defining 

'^ab,cd 



2J\f 



Sl,S2 



X 



vl/|f(^l,4). (4.7) 



n (^^^ ~ ^'=)(4 ~ ^^)(^»^ ~ w^d)(4 - ^'i) 

we propose the following expression 

^As^^ [3333] (wi, wa, W3, W4; 2^1 , • • • , 2:^^ ; 4' • • • ' 4p = 

A(0'i)[3333]({«;})vl>i2,34({«^},{^}) + 5(°'')[3333]({«;})vPi3,24({«^},{;^}) • (4.8) 

Following steps that are similar to those presented in section 13. 2| we can determine the 
coefficients in the master formula eq. ()4.8|) . The particular limits we employ are 



{II) 



T T i i 

J. i T T 

Zl, Z^,, Z[ ^ Wi,Z2 — * W2 



(4.9) 



They give 



leading to the result 



lim((T3cr3cr3a3?/^iV'iV'2V^2)^°'^^ oc {(T-icr^a^a^)^^'^^ 
lim (a-3a3a3(T3?/;iV'i?/'2V^2)^°'"^'' oc (cr|a-f 0-30-3)^°'^^ 



(4.10) 



A(°)[3333] = (w;i2W34)5x-T5(l -x)3j-f^(x) 
5(°)[3333] = (^12^34)^x^^(1 - xf^J^f\x) 
A«[3333] = -(-l)tC(wi2«;34)^x-A(l - x)^J^i^\x) 
fi(^)[3333] = -(-l)tc(w;i2W34)^x-^(l - x)ij^f ^(x) 

r(p)/ 



(4.11) 



Here, C = |y f^fjp^ and the functions J^f\x)^ with p = 0, 1 and i = 1,2 are given m 

dnHH). Note that, while we use the same notation as in section 13.21 the actual functions 
J-''f\x) and the value of C differ between the two cases. 
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4.3.2 The case n-i — ni — 2 

The wavefunction reads 



*As[ttii](wI, wl;w^^,, wir, zl, . . . , z]^^; zi, . . . , z\^^) ^ 

{a^{wl)a^{wl)ai{wi,)ai{w^)^i{zl) . . . ^i(4^)^2(4) • • -^2(4^)) 



i<j 



(4.12) 



The master formula now reads (from now on, we will label the quasi-holes consecutively, to 
avoid confusion when braiding quasi-holes; as a reminder, we will put primes on the labels 
of the spin-down quasi-holes) 

^AV^[TTii](w^I, wl;w^^,, wi,; zl, zl^^;zl, z^,) = 



with 



Sl,S2 



X 



n {zj - wDiz'^, - 

.i,j'eSi 

n {zi-wi)izi,-wi,) 

.k,i'eS2 



(4.14) 



and similarly ^'i4',23'- Note that in this case there is no natural third way to distribute the 
quasi-holes over Si, 5*2. 

To determine ^'■^'^^[TTii] and -B*-°'^^[tTii] we put — = 2 and consider the limits 



(I) 



T T J. J. 

Z^ -^2, Z^, > Z^i-, 

T TT TJ- i i i 

Z^ ^ ^1 5 ^2 ^ ^ ^3'' ^1' ^ 



(4.15) 



They give 



lim (cr-|-cr-|-crj^crj^'^i'^i'^2'V'2)^°''^^ (cr-|- cr|(T j^Uj^ 



.(0,1) 



hm (cr-f(J-f(J^cr^V'iV'iV'2V'2)^°'^^ oc ((73(73(73(73)^°'^^ . 
In the (0) channel this gives the equations 

(/) {A + 5)(°)[TTii] = (^/;i2^/;3'4')-^^^(l - x)-"^Tf\x) 



(4.16) 



(77) (v^r^A + 



1 
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5)('^)[TTii] = 



The equations in channel (1) have a similar structure. The solutions are 



(4.17) 



X 



B^%Mi] = {w,,w,,,>)-'^x'^^{l-x)-'^[^Tf\x)] 



1 — X 



(1), 



(4.18) 



Interchanging the positions of a^{w2) and a^^iw^,') gives a different basis for the four 
quasi-hole wavefunctions in this sector. With 

^AV^[TiTi](wI, w^,; wl, w],; . . . , 4, ■ ■ ■ , 4^ = 

A^'^%m\{{wm,2'M{^}. {4) + 5^'''^[TiTi]({«^})*14',32'({^}, {4) ■ (4.19) 

the coefficients are found to be 

A^%m] = (-l)^(^12'^34')-^^-^(l - X)h^f\x) - T^^f ^(^)] 

J- X 

B^%i]i] = {-l)Hwi2'W^')-'^x-^^{l- x)1[-^t'^\x)] 



-l)^C{wi2'W^i,)-ix-*{l-x)i[ri^\x) 



X 



r?{x)] 



l-x 



B^%iM] = {-i)^-c{w,2'w,,.)--^x-^^{i - x)1[-—:f\'\x)] 

-L X 



4.3.3 The case = 4 

The wavefunction takes the form 



(4.20) 



*As[tttT](wI, w^, wl, w]; 4, ■ ■ ■ , ^iv^; 4> ■ ■ ■ > ^irp = 

((7|(w|)(T|(w^)cT|(w])at(wl)^l(4) • . .V'l (4^)^2(4) • • -^2(41)) 

"^2''''^(4;4)l'^' 11(^1 -^D^IlK' -^1)^ ■ (4-21) 



X 

Prom now on, we will drop the up-arrow on the quasi-hole coordinates {wi}. Defining 



ab,cd 



2N 



Si, 52 



X 



W{zl - Wa){zl - Wi,) 

Y[{zj - Wc){zj - Wd) 
MS2 



V|>221/ T A\ 



(4.22) 
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we propose the following expression 

^{0,1) [iiii] (y^^^ W3, W4; 4, • • • , ; 4' • • • ' = 

^^°''^[TTTT](M)*i2,34({^},{4) + 5(°'^)[TTTT]({ti;})^i3,24({^},{4) . (4.23) 

To determine the coefficients in this master formula, we set A^^i = 2, A'^i = 6 and take the 
two limits 

(//) 4 ^ 4, 4 ^ 4, 4 ^ 4' 4 ^ 4' 4 ^ ^3, zl^w4 , (4.24) 

which give 

lim(a-fa|a|a|'?/'iV'i'?/'2'?/'2'«/'2V'2'?/'2'?/'2)^°'"^'' oc ((T|(jj(T|(jj)^°'"^^ 

lm(cr-f(T|Cr|Cr-fV'iV'iV'2V'2V'2V'2V'2'^2)^°'"^^ oc (cr^(J'f(J|Cr|)(°'^\ (4.25) 

Using the spin field correlators given in eq. ()B.10|) we find 

^^°^[TTTT] = (^'i2^34)^a;-*(l -a;)§^2^°)(x) 

5^'^[TTTT] = -(-l)^C(i^i2^34)^x-^(l - x)"^Ti'\x) . (4.26) 

Because the spin-up quasi-holes have to satisfy exactly the same braid properties as 
the spin-less quasi-holes for M = 0, it is not surprising at all that the functional form of the 
functions ()4.26p is exactly the same as in 1)4.111) . 

4.3.4 The case n-^ = 3, = 1 

The wavefunction is 

*As[TTTi](^I, wl, wl; wi,; z],..., zl^^;zl,, . . . , 4|) = 

{a^{wl)a^{wl)a^{wl)a^{wi)M4) ■ ■■Mzn^)M4') ■ • • V'2(4|)) 

X [¥r'\zh 4)] Ui4 - -1)^ n(4 - -i)^ ■ 

i,j i' 

X \[{wl - ^J)^ \[{w] - wi,)--^ . (4.27) 

i<j i 

In this case, we will in general have = 2r + 1 spin-up electrons and A^^i = 2r + 3 spin- 
down electrons, with r an integer. Thus, to define the functions \l/afe,cd, we will divide the 
electrons in two groups, where the first group contains (iV| — l)/2 spin-up electrons and 
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{Ni — l)/2 + 1 spin-down electrons. The second group 5*2 has the remaining [N-^ — l)/2 + 1 
spin-up electrons, and the remaining [N^ — l)/2 spin-down electrons. We can now define 



51,52 



X 



ieSi 



n (4-«^D(4 

ij'e52 



(4.2^ 



JV|+1 JV|+1 

with the normalization Af' = 2^~2^~, where the sum is over all ways of dividing the 
electrons into the two groups. Note that we again have three different ways of splitting the 
quasi-holes. However, the relation between them differs form the 'usual' relation, because 
we now have 



^23,14' = X 



12,34' 



-il-x 



WW 



13,24' 



(4.29) 



W34/ W4'2 

The master formula now reads 

"^AS^^^UKwlwlwl; wi,;zl,..., 4. • • • , 4j) = 

A^'''^[mmw})^i2Mi{^o}A^}) + B^'''\mmw})^i3M{^^ ■ (4.30) 

Specifying = 1 and A^'j^ = 3, and taking the limits 



) ^1' ' ^1) ^3' ' ^2'' ^1 



we obtain the following result 

v4(°)[TTTi] = 
5(0) [nil] = 

^(i)[llll] = 

5(i)[llll] = 



^ l5 '^3' ^ 2'' 1 



3 ' 



(4.31) 



™i:2ix-^(i-x)i^f(x) 



W34/ 



-X (1 — x)»JF-[°^(x) 



^i'4'2 



:-i)tc^^^^i^^x-A(i-x)i^i^)(^) 



W^34' 



^ 1^12^34' 



W4'2 



-x-^(l -x)ij^J^^(x) 



(4.32) 



Note that in calculating the braid relations, the 'extra' factors of W34/ and 11^4/2 are 
precisely 'compensated' by the additional factors in the relation between the \E'ab,cd''s in eq. 



4.4 Braiding relations 

The braid properties of the various four quasi-hole wavefunctions are easily evaluated using 
the transformation properties of the functions J^l^\x) as specified in section IBTsI We display 
them here for general M. 
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For braiding the neutral quasi-holes of charge 2/(4M + 3), denoted by the label '3', we 
obtain the following results 



2 iM 



KimS] (f/i2«[3333] Uu = ^ (-1)^ _^ j (4.33) 

^^[3333] ^ (f/23)?v,(g[3333] U,, = (-1)^-^ ( ^"^^^ 

*?)[3333] (f/i3)?vl/(^1[3333] U,, = (-1)*-^™ J_ "l^lj-Jj^ ) ' 
For braiding spin- full quasi-holes, of charge 1/(4M + 3), the corresponding matrices are 

= (-1)^3-3(411+3) M y U 



/ , 4 M 
U23 = (-1)^ 1 , ,,3 

^(-l)5^/r 



(-1)-^ 

l)-tv/f 



f/i3 = (-1)*-™) ( ^ J lU"^') • (4.34) 

'T — 1 sr 



These matrices are found by explicit inspection of the wavefunctions for the cases (n^ = 
3 , nj = 1) and = 4. For {n^ = 2 , rii = 2) equivalent matrices are found for situations 
where the braiding does not mix the spin-labels, such as 1 ^ 2 for ^^AslTTii] 1 ^ 3 for 

The different M-dependence between eqs. eq. (j4.33|) and (j4.34|) reflects the fact that 
for M 7^ the conformal dimension of (J3 differs from that of and cxj. For M = the 
three quasi-hole types are connected by an su{3) symmetry and the braiding properties are 
necessarily the same. 

By a simple change of basis ('gauge transformation') the matrices Uu, U23 and U13 
acquire a form which can be identified with that of the R, F and B matrices for the Fibonacci 
anyons, eq. (jl.2p . up to overall M-dependent phase factors. This establishes the suitability 
of this particular quantum register for universal topological quantum computation. 



5 Quantum group approach 

In this section, we will use the CFT-quantum group connection to calculate the braid prop- 
erties of the quasi-holes, and confirm that the results from this approach is indeed consistent 
with the results obtained here from the explicit wavefunctions for the quasi-hole states. 

The approach used here follows the lines of Slingerland and Bais, [23j, to which we 
refer for details and more references. More details on the CFT-quantum group connection 
can be found in [HHlEn]- For the su{2)k case in particular, see jSZI- More details about the 
quantum groups themselves can be found in, for instance, 

At a basic level, the connection between conformal field theory and quantum groups 
states that the braid properties of fields in the conformal field theory are the same as the 
braid properties of particles carrying a quantum group representation. Because the latter 
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are specified by the i?-matrix of tlie quantum group, one can calculate the braid properties 
of the quasi-holes in an algebraic way. In addition to the i?-matrix, one will also need to 
know the 6j-symbols, because to describe general braidings, on needs to know how to change 
between the different bases of the tensor product of three representations. This information 
is encoded in the (g-deformed) 6j-symbols, or the F-matrices, see for instance, |2I] for a nice 
review. 

To be more explicit, the F- matrices describe the basis transformation between the two 
different ways in which one can take the tensor product, or fusion, of three representations, 
or particle types. The first way is to first fuse a particle of type a with a particle of type 
b, which gives, say, a particle of type e. Finally, one fuses this particle e which the third 
particle c, with particle d as outcome. The other way of fusing particles a, b and c is to first 
fuse b and c into /, which is fused with a to give d. Pictorially, we can describe the relation 
between these two bases as 





(5.1) 



The exchange of two particles a and 6 in a definite fusion channel c is described by the 
i?-matrix. Pictorially, we have 

a\ lb a\/b 

Re'" \ =y 

\c \c. (5.2) 

In order that the F and R matrices describe consistent braiding, they have to satisfy 
consistency conditions, which go under the name of the pentagon and hexagon equations 
|22j . The F and R matrices obtained from the quantum groups automatically satisfy these 
equations. 

We can now express the braiding of quasi-holes in terms of the R and F-matrices. The 
braiding of the quasi-holes 1 and 2 (using the notation of sections El and 0]), is simply given 
by the elements of the R matrix, and depends on the fusion channel. 

The exchange of particles 2 and 3 can be done in more than one way. In our case, it 
turns out that we first need to exhange particles 1 and 2, followed by acting with the F- 
matrix, and finally, exchanging the third particle with the 'intermediate' particle. Pictorially, 
we have (for clarity, we dropped the labels a, b, etc.) 




R 




R 




This corresponds to the following form of the braid matrix U23 

(f/23)e,/ = R'/iF^^'n^R"/ 



(5.3) 



Note that we do not sum over repeated indices. Similarly, we have the following expression 
for f/i3 

(f/i3)e,/ = K^^{Ryr\F^^''rf ■ (5.4) 

We will use the expressions fl5.3|l and ()5.4|1 in the following subsections to compare the braid 
results of the previous sections with the results obtained using the quantum groups. Note 
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that by changing the relative phase of the different blocks, we can change the relative phase 
between the off diagonal elements of the matrices f/13 and f/23. Thus, only the sum of the 
off-diagonal phases is a physical quantity. Nevertheless, we reproduce the matrices exactly. 

We will not explain in detail how to obtain the F and R matrices from the quantum 
groups. For this, we refer to jSHI and a forthcoming paper jHE], in which details will be given 
on the calculation of the F and R matrices in the case of Uq{su{'i)). Here, we will be brief 
in our description, and quote the results we need from the literature. 

A quantum group associated to a Lie algebra is a g-deformation of its universal en- 
veloping algebra. At generic values of g, the representations of the quantum group are 
similar to the representations of the Lie algebra. However, if q equals specific roots of unity 
(which we will specify below in the cases of our interest), the representation theory is rather 
different. Concentration on the case Uq{su{2)) (see, for instance, [23, where the quantum 
group picture was used to calculate the braid behaviour of the quasi-holes of the Read-Rezayi 
states), we find that for q = e*+2, there are k + 1 unitary highest weight representations. 
In addition, the tensor product of such representations gets truncated in comparison to the 
tensor product of su{2). In fact, the truncated tensor products are equivalent to the fusion 
rules of su{2)k. The braid properties of the quasi-holes are in fact described by the F and 
R matrices of the quantum group at these special values of q. 

The calculation of the F matrices can be done in a similar fashion as in the case of 
ordinary groups, by first calculating the Clebsch-Gordan coefficients, and from those the 
6j-coefficients. Of course, in the quantum group case, one has to use the g-deformed raising 
and lowering operators. 

Generically, the F-matrices can be expressed in terms of the so-called q numbers \n\ , 
which are defined as 



q2 — q 2 

or, equivalently, [n\ = J2^=i <1^~^- Note that for g = 1, we have [nj = n. 

For the Read-Rezayi states with k = 3 and the paired spin-singlet states proposed by 
the authors, the corresponding value of g is g = , for which we have [Oj = 0, [Ij = [4J = 
1, [2J = [3J = and [n + 5\= - [n\ for n G N. 



5.1 The case su{2)k 

To calculate the braid properties of the quasi-holes over the Read-Rezayi states, we need 
to specify to which 5^(2)^ representation they correspond, and the possible fusion channels 
as well. Because all the braid properties are encoded by the braid properties of four quasi- 
holes and because we know the four quasi-hole correlators, we will focus on that case. The 
corresponding quantum group is Uq{su{2)), with g = The quasi-holes correspond to 

the representation I = 1 (i.e. spin i), and the two fusion channels to the representations 
/ = and / = 2. 

Let us focus on the matrix Uu first. This matrix describes the braiding of the first two 
quasi-holes, which depends on the fusion channel. For A; = 3, we have R^'^ = (— l)io and 
rI'^ = (— 1)10. This is in agreement with the matrix U12 of eq. ()3.14|) . 

11 3 2fc+l I 1 

For general k, we have Rq = — g"* = (—1)2(^+2) and R2 = g* = (— 1)2(^+2)^ which 
confirms U12 of eq. ()A.42|) . Note that the M dependence easily follows from the general form 
of the wavefunction, eq. ()A.33|1 . 
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In addition, we have for /c = 3 



r T 



Fl'' = Z \"] ■ (5-6) 



For general k, we obtain the following result 



.1,1,1 1 f -1 /[3j\ _W -I Vdt^ 



where dk = 2 cos(^^) is the quantum dimension of the fundamental representation of su{2)k. 
Now, the matrix U23 corresponds to 

iU23)eJ = R'/iF^/'rfRf , (5.8) 

with a = b = c = d= l and both e and / can take the values 0, 2. Note that we do not sum 
over repeated indices. Similarly, U13 corresponds to 

{U,s)ej = R^/iRYr'iFfy, , (5.9) 

also with a = b = c = d = l and e, / = 0, 2. 

Using these results for the R and F matrices, we easily see that the k = 3 results of 
eq. (|3.14|) and the general results of eq. ()A.42|) are reproduced exactly. 

We would like to note that the k = 3 braid matrices are, up to an overall factor, directly 
related to the braid matrices of the Fibonacci theory. Note that they are in fact related to 
the 'mirror-image' of the matrices ()1.2j) . This is related to the fact that the F matrix of the 
spin-1/2 particles in the su{2)3 theory is given by eq. ()5.6p instead of eq. ()1.2p . 

5.2 The case su{3)2 

We will now compare the results of braiding in the spin-singlet case to the results which 
can be obtained by using a quantum group picture. To do this, we will have to know the 
6j-symbols (or F-symbols) of the quantum group of su{3). The results of a direct calculation 
of these 6j-symbols for general q will be presented elsewhere ^Ei, together with the i?-matrix 
'eigenvalues' and a detailed analysis of the 6j-symbols related to cosets. Here, we will merely 
quote a small number of F-symbols and braid factors. 

For the comparison of the braid matrices with the results obtained from the quantum 
group picture, it is easiest to work with a set of parafermion fields which stay as close as 
possible to the representations used in the KZ equation, see appendix IB. II Hence, we will 
work with the representatives a = 3, 6 = 3, c = 3, and d = 8. Note that we can not take 3 as 
the last representative as well. The reason is that after fusing the first two fields, we find that 
the possible channels correspond to 6 or 3. Fusing these intermediate channels automatically 
gives us the 8. It follows that we need the following data, which can be obtained from the 
quantum group of su{3) 



where the intermediate fusion channels are 6 and 3 (in that order) and 

Rl' = (-l)A i?f = (-l)i R^ = (-1)^ i?f = (-1)^ . (5.11) 

The symbols Rg'^ and correspond to the diagonal elements of U12 as they should. In 
addition, upon using equations (j5.3p and (|5.4|) . we obtain the matrices of eq. (j4.34|) . 
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5.3 sw(3)2 parafermion correlators 



To verify the braid behaviour of the various su{3)2 parafermion correlators, we will use the 
following data 

where the intermediate fusion channels are 1 and 8 and 

Rl'^ = Rl'^ = 1 = (_l)-f = (_i)-f . (5.13) 

With this data, we exactly obtain the braid behaviour which can be derived from the cor- 
relator of four pi fields in eq. ()B.13|1 . In addition, we obtain the braid behaviour of four 
(T3 fields of eq. ()B.10|) up to an overall sign. To explain the origin of this sign, we note 
that if one expresses the SU{3)2 WZW primary in the adjoint representation in terms of the 
parafermion field a^, additional u{l) factors are needed, see eq. ()B.4|) . These u{l) factors 
give rise to the additional sign. Note that in the p sector, these u{l) factors are absent. 



A Detailed structure of the su(2)k parafermion theory 

In this appendix, we will review the structure of the parafermion conformal field theory, 
as introduced in [HI]. In fact, we will use the general setting of [SlI, so we can use the results 
to explain the details of the parafermion theory associated to S'u(3)2 in the next appendix 
as well. 

In the case of the sm(3)2 parafermions, which we will describe in the next appendix, it 
will be important to make a distinction between the chiral sectors of the parafermion theory, 
and the Virasoro primary fields. Each chiral sector contains an infinite number of Virasoro 
primary fields. However, we will often use the same notation for the chiral sector and the 
leading Virasoro primary field. For the chiral sectors we consider in this appendix, there 
will always be only one leading virasoro primary field. However, in the case of the sm(3)2 
parafermions, there will be one chiral sector containing two leading Virasoro primary fields. 

The fusion rules (which will be specified for the su{2)k parafermions below) describe 
the 'merging' of two parafermion sectors. To calculate the quasi-hole wavef unctions, we 
need to know the full details of what happens when two primary fields are brought to the 
same location inside correlators. This information is contained in the Operator Product 
Expansions (OPE's), which will be given in IA.2I in the case of the parafermions associated 
to su{2)k. 

A.l Fusion rules of the sw(2)3 parafermions 

The sectors of the su{2)k/u{l) parafermion GET (which we will sometimes denote by Z^) 
are labeled by two labels, an su{2)k label / = 0, 1, . . . , A; and the u{l) charge m, which is 
defined modulo 2k, because the u{l) theory is compactified. Thus, we write the sector (and 
the leading parafermion primary fields) as <I>^. 

The branching rules state that the only labels allowed are those which satisfy I — m = 
0mod2. In addition, we need to identify the sectors [10] $^ = It follows that there 

are ^k{k + 1) parafermion primary fields, and for each field we can choose labels with 
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/ = 0, l,...fcandme {-1 + 2,-1 + 4:,. 
of the fields are given by 

him 



I}. With labels chosen in this way, the dimensions 
1(1 + 2) 



m 

A{k + 2) ~ Jk 



(A.l) 



Note that if m is chosen outside of the range m G {—1 + 2, —l+A, ...,/}, the scaling dimension 
will be given by hi^m + where positive integer. 

Thus, for the parafermions ijji = = "^21-^? with i = 0,1, . . . , k—1 we find h^^ = 
while for the spin fields cTj = $• we find h, — 

2i G {0,1, . . . , k} have scaling dimension /i^. — -j^^. 

The fusion rules of the parafermion theory can be obtained from the fusion rules of 
the SU{2)k WZW conformal field theory. Explicitly, we have 



- 2k{k+2) ■ '^^^ 'neutral' fields Si 



$0' with 



(A.2) 



where the sum is over the range /" = |/ — /'|, |/ — /'| + 2, . . . , min(/ + V ,2k — I — V). 



Specializing to the Z3 parafermions, we will use the following standard notation 

1 = $° ^1 = $^ V^2 = $4 cri = $i a2 = ^l e = ^l. (A.3) 



The scaling dimensions are h^^ 



3' 
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and hp 



With this notation, we find the fusion rules as given in table H 



X 




0^2 


e 


^1 


^2 


o-i 


'ipi + 02 










0-2 


l + e 


^2 + 0-1 








e 


+ (Tl 


V'l + (^2 


l + e 






^1 


e 




0-2 


^2 




i'2 


(^2 


e 




1 





Table 1: Fusion rules of the parafermion and spin fields associated to the parafermion theory 
su{2),Ju{l). 

Note that the structure of the fusion rules becomes simpler if we adopt the following 
notation ifjQ = 1, ipi = ipi, ip2 = 4'2, tq = e, t\ = 02 and T2 = (J\. The labels of the fields are 
defined modulo 3. In this notation, the fusion rules are simply 

X Tj = Ti+j (A.4) 
Tj X Tj = ll>i-\-j + Tij^j 

A.2 OPE's 

In this section, we give the OPE's of the leading Virasoro primary fields^ (including the 
numerical coefficients). We will use schematic notation, writing 

A{z)B{w) = {z- w)^^-^'^-^^c%bC{w) + ... (A.5) 

^To be completely rigorous, we should view the fields as chiral vertex operators, or intertwiners. Then, 
the OPE coefficients would carry two additional labels, indicating the sectors the fields acts between. In this 
paper, we will not need this level of detail. 



22 



as 

AB = c%C (A.6) 

and restricting ourselves to the leading field in each of the fusion channels. 
We then have the following OPE's 

2 2 

^iV^l = = 1 ^2V'2 = , (A.7) 

in accordance with the results in jM] for all fields ipi and general k. The OPE's between 
parafermions and spin fields are 

^^^^ " vl'^ ^^^^ " 71^^ " 71^^ 

1 /2" 2 



while the other OPE's of the spin fields are given by, with C = ^^j p^|^p3^|^ 

0"i(Ti = -^ipi + \/2Ca2 0"icr2 = 1 + VCe: cr2cr2 = ^-^2 + V2Ccxi 
V 3 V 3 



The most interesting OPE turns out to be 



\j^il)2 + VCai a2e = ^j^ipi + VCa2 . (A.9) 



12C 

ee = 1 + Qe + J^e' , (A.IO) 

that is, the first field appearing in the 'e' channel is not the field e, but a different Virasoro 
primary field e', which has scaling dimension /i^/ = |. This result will be derived in section 
1X41 



A. 3 Spin field correlators 

In this appendix, we will explain how the four point correlators of four 'spin fields' a can be 
obtained from the four point correlators of the WZW GET as given in We will do this 
explicitly in the case of sm(2)3, and merely state the outcome in the other cases. 

The correlators which are calculated in are four point correlators of WZW primary 
fields, transforming in the fundamental representation. These fields g can be written in terms 
of the spin fields cxi and 0"2, combined with a vertex operator. Explicitly, we have gi = aie~ 
and g2 = (J2e~"f . Because correlators of the form (5'i5'25'i5'2)^^^ and {gigig2g2)^^^ are given 
in and four point correlators of the vertex operators e^"f are easily calculated, we can 
derive the explicit form of the correlators (cricr2cri(J2)^^^ and {(Tiaia2(T2)^^\ namely 

{(yi{wi)(y2{w2)cfi{w^)a2{wi))^'^'' = {wi2W34)~^x^{l - x)~^J^[°\x) 

{cri{wi)a2{w2)cri{w3)a2{w4))^^^ = (-l)^(wi2W34)"^a;^(l - x)~^CJ'f ^(x) 
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{ai{wi)cri{w2)cr2{w3)cr2{w^))^^'^ = (-1)3(^12^34) '^a; 30 (1 - x) s ^(x) 

{o-i{wi)cri{w2)a2{w3)a2{w4))^'^^ = -{-1)^ (w^Wu)'^ {I - x)^CJ^^^\x) . (A.ll) 

The constant C is related to h used in [23] via h = C"^. The phases make sure that the 
correlators have the correct behaviour when x ^ 0. In this limit, the four point correlators 
reduce to normalized two point functions. We remind the reader of the notation Wij = Wi—Wj. 
The functions J-'-^\x) can be expressed in terms of the hypergeometric functions F{a, b, c; x) 
in the following way 

^(n) , N 3 , , 1 „,1 13, ^(n) , N I7, .i„,648 , 

J- o)(x) = x--.il - x)-oFi-, --, -; x) J-f (x) = -xTo 1 - a; toF -, -; x) 

5 5 5 6 5 5 5 

1 1I37 1 1I32 

J^[^\x) =x^>(l-x)^F(-,-,-]x) y^2\x) = -2x^(1 -x)^^F(-,-,-;x) . (A.12) 

5 5 5 5 5 5 

A. 4 Further correlators 

The master formula ()3.4|1 can be used to obtain other correlation functions of parafermion 
fields, which are hard to obtain from the KZ equation, because they correspond to other 
SU(2) representations than the fundamental representations. 

These correlators are obtained by taking various limits of the master formula ()3.4|) . In 
fact, for = 3, we can obtain all the possible four point correlation functions of the spin 
fields cTi,(T2 and e, because we can fuse the spin field ai with an arbitrary number of fields ipi, 
which gives us all the possible spin fields. For arbitrary k, these methods gives us a way of 
calculating all the four point correlators of the fields of the form $^ = ^^^f,, in combination 
with an arbitrary number (which has to be allowed by fusion) of parafermion fields ipi. 

We will now use the master formula ()3.4|1 to find the correlator 
{cri{wi)a2{w2)e{w3)e{w4))''^'^\ by taking the limit {zi ^2,-22 W2,Z3 W3, w^). 
This results in 



{aiiwi)a2iw2)eiw3)eiwi)f^ = >i2"«^3/^^ (1 - a^)^ ^^^^i^) + H"' i^) (A.13) 



1 



2 



2-x 



-(0)/ X , .^(0)/ 



X 



{aiiwi)a2iw2)€iw3)e{w,)y'> = (-1)^^ ^W,2''W3/ x^ {1 - x)t [^^j-{^^(x) + J^^'>{x) 

2 11 — X 

Here, we used that C'^^ o-^ = ^J^■ We would like to note that this result is equivalent to the 
result obtained in namely 

((]•) -JL 2 142 

{cri{wi)a2{w2)e{w3)e{w4^)y = w^2^''W3^' {1 - x)~^ F{--, --, --; x) 

{ai{wi)a2iw2)e{w3)e{w^)f^ = {-ly^w^fw'^^x^l-xyl^Fi^,^,^,^) ■ (A.14) 
The equivalence of the two results follows from the fact that y/p/C = |, with p given by 

„ _ /I '^5''^ '^5''^ '^5^' (A.15) 



P(f)r2(-i)r(f)r(i; 
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and the following relations between hypergeometric functions 



2 1 4 2 . la.. 



5' 5' 5' ' 2 

a;5(l-x) 5F(-,-,— ;x) = -xio(l -x)2 
5 5 5 4 



2 -X 



1 — X 

^r(x)+^2^'^(x) (A. 16) 



.1 — X 

To obtain the OPE coefficient C^^, we expand the correlators around x = 0, with the result 

— - /'^^ X^ X^ 

Wl2'f^34 (0-l(^i^l)0-2(w2)e(w3)£(w^4))^ - 1 + — + — + o{x^) 

15 15 

7 12 17 

We find that the 1 (or ao) channel starts as x^, as was observed in As a consequence, 
we find that the OPE coefficient ^ = 0. 

Finally, we also reduced the master formula (j3.4j) to obtain the correlator 
{e{wi)e{w2)e{ws)£{w4,))^^'^\ by taking the limit {zi — > W2,Z2 — >■ ^1,2:3 — > ^3,2:4 — >■ W4). This 
results in 

{e{w^)e{w2)e{wMy^4)f^ = ^K2^34)"^a;i^(l - x)-^[(2 - x)^f (x) + (1 + x)^f (x)] 
{e{wi)e{w2)e{wMw4)f^ = ~^~^^'^, x^(l - x)-^ [(2 - x) J-f ^(x) + (1 + x)J^^'\x)] 

2{wi2W3i)5 

(A.18) 

Again, we find that the small x behaviour for the (1) or £ channel goes like xs, explicitly 

(^12^34)^ {e{wi)e{w2)e{w3)e{w4))^^'' = 1 + ^ + ^ + o(x'^) 

,4 (—1)^57 , , , , , , , , T 7x^ 261x^ / 22, 

K2W34)^ ' {siwi)eiw2)eiwMy^,)y ^ = x5 + ^ + + o(x-) (A.19) 

Using the braid properties of the functions J-'[^2^\ we find the following braid matrices 
for the correlator {e{wi)e{w2)£{w3)e{w4))^^'^^ 



4 



4 



Again, we would like to see if we can obtain the same result using the quantum group 
picture. This time, we need the R and F- matrices corresponding to / = 2 or spin-1, because 
the e particles are represented by e = $0 the coset construction. The data we need is 



-\/r -T 
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With this data, we can calculate the braid matrices using ()5.3p and ()5.4|) . with the following 
result 



-n-t 



r 



.4 / r (-r- 



We see that we almost exactly reproduced the braid behaviour of eq. ()A.20|) . The only 
difference is the additional minus signs in the off-diagonal elements of the matrices. These 
additional minus signs are a result of the fact that the conformal block in the (1) or (e) 
channel starts at a degree higher than naively expected from the fusion rules. However, these 
additional signs can be 'gauged' away, by redefining the block {eeesY^^ with an additional 
sign. Thus, we conclude that the vanishing of the coefficient ^ does show up in the braid 
matrices, but its effect can be gauged away. Note that is was to be expected that the 
additional power in the correlator would not change the braid behaviour, because there is no 
(or hardly any) freedom for change without violating the pentagon and hexagon equations. 

A. 5 Braiding relations 
From [2n] we know that 



where, in this case, we have 



^0- ^i-rfi)r^^ 



r(|)r(f) 75-1 



2 



r 



°" r(i)r(|) ci - 2r(i)r(|) ^^'^'^ 

In addition, we have the following expression for C, which is related to h used in [23] 
For general fc, these relations become 



cl = -cl 



In addition, we have the following expression for h, see 



^2 _ C*? _ 1 r(fc_^2)^(fc+2)^ (fc+2) Cl _ I _ I ^0 „i . . . 
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The transformation behaviour of the JF/^Vx) under x t— > -r^ is as follows 

= _(_i)-2(^(i_a;)2(iT2)J-(°)(x) 

1 — X 

•^i'^(lT^) = -(-1)^(1 -a;)^[J-f^(a;) + 

jr(i)(^^) = (_i)2(feW(l - J-o^'Vx) . (A.28) 

1 — X 

For X I— > i we have 

X 

^f^(^) = (-l)W2)a;5a?W[C°^f)(x) + ^(x)] (A.29) 

^2^°^^) = -(-l)5{^xW2)[C°(j-f^(x) + J-f^(x)) + (-l)^C°(^f^(x) + J-f^(x))] 
J-f ^i) = -{-ir^x^[C'oJ^i^\x) + (-l)^C!^|'^(x)] 

^2^'^-^) = (-l)-^xW2)[Col(^(°)(x) (X)) + (-l)^Cl(j-J'^(x) + J-2^'^(X))] . 

A. 6 More general k results 

From [22], we can extract the following correlators for general k 

(ffi(wi)fTfc_i(w2)ai(w3)afc_i(u;4))^°^ = {w^Ws^y^ x^ {1 - x)-^^f ^(x) (A.30) 
(ai(wi)afc_i(w2)ai(w3)afe_iK))^^^ = i~l)^^iwuW34r^x^{l ~ x)-^CJ^i^\x) 

{cTl{Wi)cTi{w2)cTk-l{w3)cTk^l{w4))^^'^ = (- 1) ^ (W12W34) " (1 _ 2X J^W 

{cxi{wi)cxi{w2)cxk^i{w3)ak-i{wi)f^ = -{-l)T^{wi2Wur^x^^{l-x)^CJ^^^\x) , 
where the J-'-^'^ are now k dependent 

J^f^fx) = x"^(^(l - x)^(^F(- , , :x) 

W a; = a;2(fc+2) 1 OTTi?^ - x) 

J^\^\x) = x^o^(l - x)'^o^F(-^—, x) 
^ ^ ' ^ ' ^k + 2' k + 2' k + 2' ' 

J^^^\x) = -2x^(^(1 - x)^(^F(-^,-^,-^;x) . (A.31) 
^ ^ ' ^ ' ^k + 2' k + 2' k + 2' ' ^ ' 



From this, we find the following OPE coefficients 

1 

7^ 
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The wavefunction of four quasi-holes and N — {r + l)k — 2 electrons (with r a positive 
integer) can be written in the following form 

^(°'^^(wi,W2,W3,W4;zi, ...,zn) = {(Jiiwi)ai{w2)ai{ws)ai{w4)tlJi{zi) ■ ■ ■ ^i(^iv))^°'^^ x 

. . 3 2 1 -w—r , , 1 -I — r fcM+2 , . , 

X (wi2t(;34)M'=M+2)3; KkM+1) (1 _ 2;)fe(fcM+2) YYwi - Zj)k ^^{Zi - Zj)~^ . (A.33) 

Again, we have the following master formula 

*(°'^)(Wl,W2,W3, W4;^l, ...,Zn) = 

^(°'^n{^})*(12)(34)({^},{4) + ^^°''H{^})*(13)(24)(M,{4)- (A.34) 

To specify the functions ^(i2)(34) and ^(i3)(24), we divide the electrons into k groups, namely 
Si^{i,i + k,...,N-{k-2)+i} for i = 1, . . . , and 5"^ = {j,j + k,...,N-{2k-2)+j} 
hr j — k — 1, k. Setting M = for simplicity, we have 

1 



*(i2)(34) ^jfYii n ^^:>~ ~ "^2) n (^i' ~ - n 



{Si} ~jeSk-i j'eSk i=i 



^(13)(24) 



1 



jGSk-1 j'&Sk i=l 



where the sum is over all in-equivalent ways of dividing the electrons in k groups and A'' 
We will consider the following two hmits in the case oi N — 2k — 2 electrons 



Zi - 


Zl 


I — 


Zl - 


W2 




^3 - 


Zk 


3 = 


Zk - 


-> W4 





Zi ^ Zl i — 2, . . . , k — 1 

Zl W3 

Zj ^ Zk j = k + l,...,2k-2 

Zk^W4,. 

(A.36) 

On the one hand, the master formula reduces to a form containing the following two cor- 
relators (cricrfc_icricr/;_i) and (cricri(Tfc_i(jfc_i) in the limits (/) and (//) respectively. On the 
other hand, we find 



((^-1) 



^2 



lim^(i2)(34) = W14W32W42 ^ lini ^(i2)(34) = (A.37) 



((^-1)!)' 2.-2 



lim*(i3)(24) = lini ^'(i3)(24) = TfeTI «^14W32'U;34 

(/) {II) K 

The functions A^'^ and B^'^ again follow 

^(0) = {wi2WuTx^{l - Xyj^f'\x) 
S(0) = -{W12W34TX^-\1 - xy+^J^i''\x) 
^(1) = -(-l)42C(wi2W34)"x''(l - Xyj^l^\x) 
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5(1) = {-l)T^^C{wuW34Tx'^'\l - xY+^J^^^\x) , (A.38) 
where we introduced the following notation 

_ 2k + l 3M f^- ^ _ M 

^~ 2{k + 2)~ 2{kM + 2) ^~ 2{k + 2)^ 2{kM + 2) ^ ~ ~2{kM + 2) ' ^ ' 

In this derivation, we used that the OPE coefficients for the parafermion fields ipi are 
given by [21] 

2 r(/ + r + i)r(fc-/ + i)r(A:-r + i) 



from which it follows that 

1=1 

In addition, we find the following braid matrices 



2fc+l 

M I 1—1 ^ 

= (-1) 2{fcM+2) 



l)2(fc+2) 

1 

(^_l)2(fc+2) 
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2(^+2) 2{fcM+2) / X (-l)"fc+2^d2 _ I 

4 W-l)^v^d2r^ _i 



Qifc y-Vfifc-l (-l)fc+2 



where 4 = (C^)"- 



B Detailed structure of the su(3)2 parafermion theory 

We now turn the su{3)2 parafermions, as introduced by Gepner in [25 • This theory arises 
upon factoring two free fields from an SU{3)2 WZW theory, namely, it is the sti(3)2/[u(l)]^ 
coset CFT. In this appendix we present fusion rules, OPE's, 4-point correlation functions 
and braiding properties. 



B.l Fusion rules 

The su{3)2 parafermion theory has 8 chiral sectors, which we label as 

{l,V'l,^2,V^12,Crt,f^i,CT3,P} • 

It is important to realize that each of these sectors contains an infinite number of Virasoro 
primary fields. We sometimes use the same notation for the parafermion sector (say V'l) and 
for the leading Virasoro primary field (ipilz)). We shall see that in the sector denoted as 
p there are two independent leading Virasoro primaries, which we shall write as Ps{z) and 

Pc{z). 

At the level of parafermion sectors, the merging of two sectors is expressed via fusion 
rules, which we present in this section. Our computations in section HI require more detailed 
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information at the level of (primary) fields. The latter is contained in the OPE's that we 
present in the next section of this appendix. 

The fusion rules of the eight parafermion sectors can be derived from the coset descrip- 
tion of this theory. For our purposes, we will use the defining coset sm(3)2/['u(1)]^, as given in 
[SI]. The parafermion sectors are labeled by an su{3) label A = (Ai, A2) (when labeling the 
representations in terms of the dimensions, we would have 1 = (0, 0), 3 = (1,0), 3 = (0, 1), 
6 = (2, 0), 6 = (0, 2) and 8 = (1, 1)) and two u{l) labels A = (Ai, A2). 

There are various restrictions on the labels. First of all, we have the branching condition 
Ai + 2A2 = (Ai + 2A2) mod 3. The label A is only defined up to 2 times (in general, k times) 
the root lattice of su{3), which means the following sectors are identified 



A 

Ai+4,A2-2 



Ai,A2 



Ai-2,A2 



(B.l) 



In addition, there are other identifications, which follow from the structure of the affine Lie 
algebra s-u(3)2, see 



(Ai,A2) 
(Ai,A2) 



^(2-Ai-A2,Ai) ^ ^(A2,2-Ai^ 



■A2) 



(Ai+2,A2) 



■ (Ai,A2+2) 



(B.2) 



From these rules, it follows that there are indeed eight different sectors, or 'parafermion 
fields', as mentioned above. The fusion rules follow from the general rule 



A"eAxA' 



A+A' 



(B.3) 



The fusion rules can now easily be derived from the following set of 'representations' 
of the eight parafermion sectors, which close under fusion 



$(0,0) 
(0,0) 



(1,1) 
(-1,2) 



$(0,0) 
(2,-1) 



^2 



$(0,0) 
(1,-2) 
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$(0,0) 
(1,1) 



(1,1) 
(2,-1) 



CT3 = $ 



(1,1) 
(1,1) 



(1,1) 
(0,0) 



(B.4) 



Using the only non-trivial fusion rule of the su{3)2 fields (0, 0) and (1, 1) (corresponding to the 
one and eight dimensional representation, respectively), namely, (1, 1) x (1, 1) = (0, 0) + (l, 1) 
and the identifications ()B.1|) . we find the fusion rules as given in table |21 



X 


0"! 




f^3 


P 




i^2 


^12 


cr 1 


1 + P 
















1pl2 + CTs 


1+P 












0-3 




1p2 + 0"T 


1 + P 










P 


1p2 + C^t 


^1 + 0-i 


^12 + 0"3 


1+P 








^1 




p 






1 






^"2 


P 


0-3 






i^i2 


1 




-012 






P 


(^3 


^2 




1 



Table 2: Fusion rules of the parafermion and spin sectors associated to the parafermion 
theory sm(3)2/[m(1)]^ introduced by Gepner [Slj. 

Of course, these fusion rules can also be derived from the S'-matrix and the Verlinde 
formula. 
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B.2 OPE's 

To our knowledge, the OPE's of the leading Virasoro primary fields in the various sectors of 
the su{3)2 parafermion theory have not been presented in the literature. 

We have determined the leading terms in the OPE's of these Virasoro primaries. We 
observed a Z3 symmetry relating (cxj^, (T|, (T3) and {ipi,ip2, V^ipu)- To streamline notations, 
we therefore write ai = ai, a2 = cr^, and ip^ = \/2'ipi2. We also employ Virasoro primaries pi 
and Pi, i = 1,2,3. These fields are all linear combinations (specified below) of the fields ps 
and pc- The scaling dimensions of the leading fields are 

113 

K = 2^ ^'^ " To ' ^'^ " 5 ■ ^^'^^ 

For fixing the details of the OPE's, especially those involving the fields p^, Pc, we 
have relied on various contractions of the master formulas such as eq. ()4.8|1 . ()4.13|1 . etc., 
presented in section 14.31 These formulas provide expressions for correlators of four spin 
fields (Tj plus an arbitrary number of parafermions tpj. By fusing some of the ai with the ipj, 
one produces various combinations of the fields c ; in the end this gives enough information 
to uniquely fix the OPE's. [We remark that the logical structure of our reasoning has been 
quite delicate: we have relied on the general 'master formula' structure of correlators and 
on the 'seed' provided by the explicit correlators provided by the KZ paper and by some 
of the simplest OPE's, set by the parafermion fusion rules. Combined these turn out to be 
strong enough to fix both the coefficients in the master formulas and the OPE's of all fields 
involved.] In section lB.41 below we explicitly mention some of the contractions we used and 
we provide some additional correlation functions of the parafermion theory. 

Employing the same schematic notation as in section eq. ()A.5|) . ()A.6|) . we have 
the following OPE's 



1 for i = j 
1 

for t ^ j, t ^ k, j ^ k 
Pi for i = j 

-^cTfc for i ^ i ^ k, j ^ k 



1 + V2Cpi for i = j 
1 

7! 



+ V-3C o-fc for i ^ j, i ^ k, j ^ k 



(jj for i = j 

ipj + . . . for i = i 

+ ■ ■ ■ ioii^j 



2Co'j + . . . for i = j 



(B.6) 
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Note that we have written separate equations for fusing fields from the p and a sectors into 
the (0) channel or the (1) channel. 

The pi and pi can all be expressed in two independent Virasoro primaries ps and pc 
according to 

Pi = ^(Pc - VSp^) , P2 = ^(Pc + VSps) , p3 = -pc, 

Pi = ^(v^Pc + Ps) , P2 = ^(-V^Pc + ps) , P3 = -ps , (B.7) 

The list of basic OPE's is then completed by 

PsPs = 1 - V2Cpc , PsPc = -V2Cps , PcPc = 1 + V2Cpc (B.8) 

giving 



PiPj = 1 - V 2Cpi for i = j 

= -il - V2Cpk for i ^ j, i ^ k, j ^ k 



PiPj = 1 + V 2Cpi for i = j 

= ~l + V2Cpk for J, t^k, j^k (B.9) 

In all these expressions, the constant C has the value given below eq. ()4.1H] . 
B.3 Spin field correlators 

We now present the four-point correlation functions for the spin fields in the parafermion 
theory. They are obtained by factoring the four point functions of fundamental fields in the 
SU {3)2 WZW model, as given by by factors associated to the spin and charge bosons 

{a^iwi)a^iw2)a^iw3)aiiw,))^'^ = iwi2W3,)~h^ (l - ^(x) (B.IO) 

{a^{wi)a^{w2)cri{w3)a 1(104))^^'^ = -(-l)tc(wi2ti'34)"^x^ (1 - x)'^J^i^\x) 

(crt(wi)cr|(w2)at(w3)cr|(w4))^°^ = (-1)^ (wuWa^)'^ {1 - x)^i J='^°\x) 

{ai{wi)ai{w2)(T^{w3)ai{w4:))''^^ = (-l)^C(wi2W34)~^x^(l - x)^i J='i^\x) 

{a3{Wi)a3{w2)cT3{w3)a3{w4))^°'^ = {Wi2W34)~^xi^{l - x)'^ [J^[°\x) + J^2°^(x)] 

{cr3Ma3{w2)cr3{w3)cr3{w,))^'^ = -{-l)iC{wuW34)^'^x^o {1 - x)^J'i'\x) + J'^'\x)] , 
where x is as usual x = ("'i~"'2)(w3-w4) r^j^^ functions J^^^\x) are now given by 

Tf\x) = - a;)i^F(l, -i, | x) Tf\x) = 1x^(1 - ^, ^; x) 

T^\x) =x^^{l- % \- x) t!^\x) = -3xii(l - x)^^F{\ % % x) . (B.ll) 

5 5 5 5 5 5 
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B.4 Further correlators 



The master formulas developed in section 1^731 can be used to gain insight into the OPE's and 
correlation functions involving the fields in the p sector of the sm(3)2 parafermion theory. 
In section |R2] we already displayed some of the OPE's satisfied by ps and pc- In deriving 
these we proceeded as follows. Having defined the combinations pi and pi through the OPE's 
((TjCTj)*^^^ = \/2Cpi and (■^/^jCTi)^^^ = pi, we observed that a specific contraction of the master 
formula eq. ()4.13|) shows that {(Jiaipi) vanishes. This implies that pi and pi are orthogonal 
(have a vanishing 2-point function). Exploiting the symmetry among the i = 1,2,3 labels 
leads to the parametrization eq. ()B.7|) . All this leaves some freedom in the self-couplings of 
the Pc and ps. To fix these, we observed that yet another contraction of the master formula 
implies that the three point function {P1P2P3) vanishes. This contraction arises from the 
formula 



((Ji ( U7i ) 0-2 ( ) 0-3 ( 3 ) V^l ( ) ( ) ^3 ( ^3 ) ) 



X 



{Zi - W2){Z2 - W3){Z3 - Wi) + {Zi - W^){Z2 - Wi){z^ - W2) 
[{Zi - W2){Z2 - W3){Z3 - Wi){zi - W3){z2 - Wi){z2, - 102)]^ 



(B.12) 



in the limit where Zi Wi for i = 1, 2, 3. 

The vanishing of {P1P2P3) implies that the combination of pc and ps featuring in the 
(1) channel of the fusion product of pi and p2 is orthogonal to ps and thereby proportional 
to P3. A final contraction yielding {a2(J2pipi) is then used to fix the normalization, giving 
eq. dHSl). 

Correlation functions involving one or more fields in the p sector are easily generated as 
suitable contractions of the various master formulas. Examples are the following four point 
functions 



(0), 



{Pi{wi)pi{w2)pi{w3)Pi{wa))^^^ = 

(pi(Wl)pi(w2)pi(u^3)Pl(w4))^^^ = 

-{-l)IC{wi2Wu)'^^x^^{l - x)-l{l -x + x^) [j^i^\x) + J^^^\x) 

(pi(Wl)pi(w2)p2 (1^3)^2(^4))^°^ = 



-(W12W34) 5X15(1 -a;) 3 

{pl{Wi)pi{w2)p2{m)p2{WA)Y^^ = 
3 C 6 8 

-(-1)5 j(Wl2W34)~^xT5(l - X 



[i-xy 



J^i'\x) + 



X 



:i-x) 



2:Ff\x) 



(B.13) 



1 + 



{i-xy 



J'\'\x) + 



X 



1-x) 
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B.5 Braiding relations 



We list the transformation properties of the functions T{2 specified in eq. flB.ll|) 

under transformations (i) x — > 1 — x, (ii) x t^-, (iii) x ^ -. 



For a; —> 1 — a; we have [22] 



T'^l{\-x) = yClT^l{x) (B.14) 



with 



°" r(f)r(t)' a - 3r(i)r(f)- ^^'^'^ 



We observe that 



g^ = gi7Co = Q r|||r3(3) > = Co^c = -v^ . (B.ie) 



For X — "T"^ we have 

1 — X 



Finally for x 



Tf\^) = {-l)-Hl-x)\Tf\x) 
\ — X 

•^f (t^) = - x)^\-xTf\x) + (1 - (x)] 

Ti'\-p^) = i-l)^^{l-xf^T['\x) 

1 — X 

^i'^(T^) = (-1)^(1 -x)t[-x^W(x) + (1 -x)^(^)(x)]. (B.17) 



(1) = (-l)-*x^[C°^f (x) - (-l)tC7?^«(x)] 
^W(i) = (-l)-^xi[-(-l)tCo^^f (x) - {-l)-y'Cl4'\^))] 

dj 

= (-l)-Ax^[C° (x^f^(x) - (1 -x)^f (x)) 



X 



-l)lcl{xT['\x)-{l-x)T^^\x))] 



= (-l)-Axi[-(-l)tCo^ (xJ-f)(x) - (1 - x)4'\x)) 

X 



i-l)-^Cl {xJ^i'\x) - (1 - x)4'\x))] . (B.18) 
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